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Abstract 

The characteristics of the memory of accelerated motion in Min- 
kowski spacetime are discussed within the framework of the nonlocal 
theory of accelerated observers. Two types of memory are distin- 
guished: kinetic and dynamic. We show that only kinetic memory is 
acceptable, since dynamic memory leads to divergences for nonuniform 
accelerated motion. 
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1 Introduction 



The special theory of relativity is based on two basic postulates: Lorentz 
invariance and the hypothesis of locality Lorentz invariance refers to a fun- 
damental symmetry principle, namely, the invariance of basic physical laws 
under inhomogeneous Lorentz transformations. In practice these laws of na- 
ture involve physical quantities measured by inertial observers in Minkowski 
spacetime. An inertial observer always moves uniformly and refers its ob- 
servations to the fixed spatial axes of an inertial frame; it can be depicted 
by a straight line in the Minkowski diagram and represents an ideal; in fact, 
physical observers are all effectively accelerated. For instance, one can imag- 
ine the influence of radiation pressure on the path of a cosmic particle. In 
general, the acceleration of an observer consists of the translational accel- 
eration of its path as well as the rotation of its spatial frame. Observers 
with translational acceleration are therefore represented by curved lines in 
the Minkowski diagram. As an example of a rotating observer, consider a 
uniformly moving observer that refers its observations to spatial axes that 
rotate with respect to the spatial frame of the underlying inertial coordinate 
system. The hypothesis of locality refers to the measurements of realistic (i.e. 
accelerated) observers: such an observer is postulated to be equivalent, at 
each event along its worldline, to a momentarily comoving inertial observer. 
The origin of this assumption can be traced back to the work of Lorentz in 
the context of his classical electron theory p] ; later, it was simply adopted 
as a general rule in relativity theory ||. 

Along its worldline, the accelerated observer passes through a continu- 
ous infinity of hypothetical momentarily comoving inertial observers. Stated 
mathematically, the translationally accelerated observer's curved worldline 
is the envelope of the straight worldlines of this class of hypothetical iner- 
tial observers. Therefore, the hypothesis of locality has two components: (i) 
the assumption that the measurements of the accelerated observer must be 
somehow connected to the measurements of the hypothetical class of mo- 
mentarily comoving inertial observers along its worldline and (ii) that this 
connection is postulated to be the pointwise equivalence of the accelerated 
observer and the momentarily comoving inertial observer. The latter means 
that the acceleration of the observer does not directly affect the result of 
its measurement; devices that obey this rule are called "standard". Thus 
the hypothesis of locality is a simple generalization of the assumption that 
the rods and clocks of special relativity theory are not directly affected by 
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acceleration 0. 

What is the physical basis for the hypothesis of locality? It is difficult to 
argue with part (i) of this hypothesis, since the fundamental laws of (non- 
gravitational) physics have been formulated with respect to inertial observers 
and hence the measurements of accelerated observers should be in some way 
related to those of inertial observers. On the other hand, part (ii) can only be 
valid if the measurement process occurs instantaneously and in a pointwise 
manner. That is, (ii) is appropriate for phenomena involving coincidences of 
classical point particles and null rays. Classical waves, on the other hand, are 
extended in time and space with a characteristic wave period T and a corre- 
sponding wavelength A, respectively. Imagine, for example, the measurement 
of the frequency of an incident electromagnetic wave by an accelerated ob- 
server; at least a few periods of the wave must be received by the observer 
before an adequate determination of the frequency would become possible. 
Thus this measurement process is nonlocal and extends over the worldline 
of the observer. The observer's acceleration can be characterized by certain 
acceleration lengths L given by c 2 /g and c/Q for translational acceleration 
g and rotational frequency Q, respectively. The nonlocality of the external 
radiation is thus expected to couple with the intrinsic scales associated with 
the acceleration of the observer. 

Classical wave phenomena are expected to violate the hypothesis of local- 
ity. The scale of such violation would be given by X/L = T/(L/c), where L/c 
is the acceleration time. The hypothesis of locality will hold if A is so small 
that the incident radiation behaves like a ray, i.e. in the eikonal (or JWKB) 
limit such that X/L —>■ 0; alternatively, L can be so large that deviations of 
the form X/L would be below the sensitivity threshold of the detectors avail- 
able at present. Consider, e.g., laboratory experiments on the Earth; typical 
acceleration lengths would be c 2 /g (B ~ 1 lyr and c/VL® ~ 28 AU, so that for 
essentially all practical purposes one can ignore any possible deviations from 
locality at the present time. In this way, we can account for the fact that 
the standard theory of relativity is in agreement with all observational data 
available at present. As a matter of principle, however, it is necessary to 
contemplate generalizations of the hypothesis of locality in order to take due 
account of intrinsic wave phenomena for realistic (accelerated) observers. 

All of our considerations in this paper are within the framework of classi- 
cal field theory; nevertheless, it is necessary to remark that quantum theory 
is based on the notion of wave-particle duality, and so an adequate treatment 
of classical wave phenomena is a necessary prelude to a satisfactory quantum 
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theory. 

To proceed, we consider the most general extension of the hypothesis of 
locality that is consistent with causality and the superposition principle. A 
nonlocal Lorentz-invariant theory of accelerated observers has been developed 
along these lines || f|, [|, |H and is presented in Section 2. The theory involves 
a kernel that depends primarily on the acceleration of the observer; that is, 
the measurements of the observer depend on its past history of acceleration. 
The main physical principle that is employed in the nonlocal theory for the 
determination of the kernel is the assumption that an intrinsic radiation field 
can never stand completely still with respect to an accelerated observer; this 
statement involves a simple generalization of a property of inertial observers 
to all observers. Thus the accelerated observer is endowed with memory, 
and the past affects the present through an averaging process, where the 
weight function is proportional to the kernel K(t,t'). It turns out that 
the kernel K cannot be completely determined by the theory presented in 
Section 2. An additional simplifying assumption is therefore introduced in 
Section 3: K(t, t') must be a function of a single variable. Two cases are then 
considered: (1) K(t,t') = k (r') and (2) K(t,t') = k(r-r'). We show that 
case (1) — i.e. the kinetic memory case — has acceptable properties that are 
described in Section 3. Case (2), i.e. the dynamic memory case, is treated in 
detail in Sections 3 and 4, where it is shown that the kernel function k can be 
unbounded even if the observer's past history has constant velocity except 
for one episode of smooth translational acceleration with finite duration. 
Specifically, we study the measurement of electromagnetic radiation fields by 
an observer that undergoes translational or rotational acceleration that lasts 
for only a finite interval of its proper time. After the acceleration is turned 
off, the observer measures in addition to the regular field a residual field that 
contains the memory of its past acceleration. This leftover piece is a finite 
constant field (kinetic memory) in case (1); however, it is time dependent 
(dynamic memory) in case (2). We rule out the latter case, since we prove 
that the measured field could diverge under certain reasonable circumstances. 
We are thus left with a unique theory that involves kinetic memory. An 
important aspect of our nonlocal ansatz is that the kernel induced by the 
acceleration of the observer depends on the spin of the radiation field under 
consideration. In particular, the kernel vanishes for an intrinsic scalar field, 
i.e. such a field is always local. As discussed in Section 5, our theory therefore 
rules out the possibility that a pure scalar (or pseudoscalar) field exists in 
nature. This conclusion is in agreement with available experimental data. 



4 



The nonlocal theory therefore predicts that any scalar particle would have to 
be a composite. Section 5 contains a brief discussion and our conclusions. A 
detailed discussion of the observational consequences of the nonlocal theory 
is beyond the scope of this work. In the following, we use units such that 
c—1, i.e. the speed of light in vacuum is unity. 

2 Accelerated observers and nonlocality 

The measurement of length by accelerated observers involves subtle issues 
in relativity theory that have been investigated in detail ]7| || |J; for our 
present purpose, the main result of such studies is that an accelerated frame 
of reference, i.e. an extended coordinate system set up in the neighborhood of 
an accelerated observer, is of rather limited theoretical significance. We shall 
therefore refer all measurements to an inertial reference frame in Minkowski 
spacetime. 

Imagine a global inertial frame with coordinates x = (t, x) and the stan- 
dard class of static inertial observers with their orthonormal tetrad frame 
A^ = 5£, where Aj^ is the temporal direction at each event and A^, 
i = 1,2,3, are the spatial directions. The hypothesis of locality implies that 
an accelerated observer is also endowed with a tetrad frame A^Jr), where r 

is the proper time along its worldline. For each r, A(* x(r) coincides with the 
constant tetrad frame (related to A/ \ by a Lorentz transformation) of the 
momentarily comoving inertial observer. We note that dX^/dr = (pj 3 ^^, 
where 4> a/ 3 = —4>/3 a is a tensor such that 4> 0i = (g)j and 0^ = e^kiStyk- 
Here g(r) is the translational acceleration of the observer and f2(r) is the 
rotational frequency of its spatial frame. Each element of the acceleration 
tensor <fi a p is a scalar under the inhomogeneous Lorentz transformations of 
the background spacetime. We assume throughout that the acceleration is 
turned on at r = r and will in general be turned off at T\ > r . 

Let f^ u represent an electromagnetic radiation field as measured by the 
standard set of static inertial observers. According to the hypothesis of 
locality f a p = //^A^A^, i.e. the projection of the field on the instantaneous 
tetrad frame, would be the field measured by the accelerated observer. On 
the other hand, let F a p(r) be the true result of such a measurement. Taking 
causality into account, the most general linear relationship between F a p{r) 
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and f a p{r) is 

F a0 (r) = f a p{r) + I K a(3jS (r, r')P\r') dr> . (1) 

This relation refers to quantities that are all scalars under the Poincare group 
of spacetime transformations of the underlying inertial coordinate system. 
We note that the magnitude of the nonlocal part of equation (JT|) is of the 
form X/L if the kernel is proportional to the acceleration of the observer. 
It follows from Volterra's theorem that in the space of continuous functions 



the relationship between F and / is unique [10, 11]; this theorem has been 
extended to the Hilbert space of square-integrable functions by Tricomi (12 . 

The basic ansatz ([]]) is consistent with an observation originally put for- 
ward by Bohr and Rosenfeld that the electromagnetic field cannot be mea- 
sured at a spacetime point] in fact, an averaging process is necessary over 



a spacetime neighborhood [13, 14]. In the case of measurements by inertial 



observers envisaged by Bohr and Rosenfeld |TB|, 13], there is no intrinsic tem- 
poral or spatial scale associated with the inertial observers; therefore, one 
can effectively pass to the limiting case of a point with no difficulty as the 
dimensions of the spacetime neighborhood can be shrunk to zero without any 
obstruction. For an accelerated observer, however, the intrinsic acceleration 
time and length need to be properly taken into account. Hence the nonlocal 
ansatz ([[]) may be interpreted in terms of a certain averaging process over 
the past worldline of the accelerated observer. 

To determine the kernel K, let us first mention a basic consequence of 
the hypothesis of locality for a radiation field. Imagine plane monochro- 
matic electromagnetic waves of frequency u propagating along the z-axis 
and an observer rotating uniformly about this axis with frequency f2 in the 
(x, y)-plane on a circle of radius p in the underlying inertial reference frame. 
We find from f a p = f^X^X^ that according to the rotating observer the 
frequency of the wave is Cj = ^{uj =f Qq), where 7 is the Lorentz factor corre- 
sponding to the speed pQo of the observer and the upper (lower) sign refers 
to incident positive (negative) helicity radiation. This result has a simple 
intuitive interpretation: In an incident positive (negative) helicity wave the 
electric and magnetic field vectors rotate with frequency uo{— uo) about the 
direction of propagation of the wave. As seen by the rotating observer, the 
field vectors rotate with frequency u — Q (—u — Q ) with respect to the 
inertial temporal coordinate t; moreover, the Lorentz factor simply accounts 
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for time dilation dt = ^dr. It follows that a positive helicity incident wave 
can stand completely still with respect to all observers rotating uniformly 
with frequency Qq = uj. In terms of energy, we have E = ^(E—cr- f2o), 
where cr is the spin of the incident photon. More generally, for oblique inci- 
dence E = ^{E — HMQo), where M is the multipole parameter such that hM 
is the component of the total (orbital plus spin) angular momentum along 
the z-axis. This is an example of the general phenomenon of spin-rotation 
coupling; various aspects of this effect and the available observational evi- 
dence are discussed in |15|, [Tj| [17|, [18|, [19|, |(J. Again, the incident wave can 



theoretically stand completely still for all observers rotating with frequency 
Qo such that uj = MQ . Let us recall here a fundamental consequence of 
Lorentz invariance, namely that a radiation field can never stand completely 
still with respect to an inertial observer. That is, an inertial observer can 
move along the direction of propagation of a wave so fast that the frequency 
uj = ju(l — 0) can approach zero but the mathematical limit of uj = is 
never physically achieved, since the observer's speed cannot reach the speed 
of light in vacuum (/3 < 1). Therefore, for an inertial observer uj = implies 
that uj = 0. On the other hand, while we find that the hypothesis of locality 
predicts that a circularly polarized wave can stand completely still with re- 
spect to a uniformly rotating observer, this possibility can be avoided in the 
nonlocal theory by an appropriate choice of the kernel. 

To implement the requirement that a radiation field can never stand com- 
pletely still with respect to any observer, we assume that if F a p(r) turns out 
to be constant in equation (|l|), then must have been originally constant 
just as in the case of inertial observers in the standard theory of relativity. 
It is convenient to replace the tensor f^ u by a six-vector /, with electric and 
magnetic fields as components, and introduce the "Lorentz" matrix A such 
that / = Af. Then for constant fields / and F, equation (TO) takes the form 

F = A(r)/ + f K(r, r')A(r')/ dr', (2) 

where for r = ro, the matrix Ao := A(tq) is constant and F = Ao/. Thus 
in the nonlocal theory the kernel K should be determined from the Volterra 
integral equation 

A = A(r)+ f K(T,r')A(r')dr'. (3) 

J TO 
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It follows from Volterra's theory (see Appendix A) that to every kernel K 
corresponds a unique resolvent kernel R(t, t') such that 

A(r) = A + f R(r,r')A dr'. (4) 

J TO 

Therefore, only the integral of the resolvent kernel is completely determined 
by our physical requirement 

f R(r,r')dr' = A(r)A 1 - J, (5) 

J TO 

where / is the unit matrix. It is clear at this point that given A(r), relations 
(0)-(|§) are not sufficient to determine the kernel K uniquely To proceed 
further, other simplifying restrictions are necessary on K or R, 



f(r) 



F(r)+ f R(r,T')F(T')dr'. (6) 



This must be done in such a way as to preserve time translation invariance 
in the underlying inertial coordinate system. 

Let us finally remark that for a scalar field, A(r) = 1 and equations (^)- 
@ simply reduce to the requirement that K(t, t') must have a vanishing 
integral over r' : r — > r. That is, the connection between the kernel and 
the acceleration of the observer disappears. This circumstance is further 
discussed in Section 5. 



3 Memory 

It is necessary to introduce simplifying assumptions in order to find a unique 
kernel K. We therefore tentatively postulate that K is a function of a single 
variable. There are two reasonable possibilities: 

K{t,t') = k {r') (easel) 

and 

K{t,t') = k{r -r'); (case 2) 

in either case, the basic requirement of time translation invariance in the 
background global inertial frame is satisfied. 



8 



3.1 Kinetic memory 

In case (1), the kernel ko corresponds to a simple weight function that can 
be determined by differentiating equation @, 

The kernel ko is thus directly proportional to the acceleration of the observer. 
A significant feature of this kernel is that once the acceleration is turned off 
at r = Ti, then for r > r%, 

F{r) = f{r)+ I" k (r')f(T') dr'. (8) 

J TO 

There is therefore a constant memory of past acceleration and the field F 
satisfies the standard field equations in the inertial frame. That is, the field 
equations are linear differential equations and the addition of a constant 
solution is always permissible but subject to boundary conditions. In terms 
of actual laboratory devices that have experienced accelerations in the past, 
such constant fields as in equation (|8p would be canceled once the devices 
are reset. Thus case (1) involves simple "nonpersistent" memory of past 
acceleration; therefore, we call ko the kinetic memory kernel. 

It is interesting to note that our basic integral equation (H) together with 
the kinetic memory kernel (0) and an integration by parts takes the form 

F(r) = F(t ) + I Adf, 

J [tq,t] 

so that dF = Adf along the worldline of the accelerated observer. 

3.2 Dynamic memory 

The second case involves a convolution type kernel K = k(r — r'). It follows 
(see Appendix A) that in this case the resolvent kernel is of convolution type 
as well, R = r(r — r'). Thus equation (||) can be written, after expressing 
the left side as the area under the graph of the function r from the origin to 

r — r = t, as 
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The kernel k is then given by (cf. Appendix A) 



k(t) = -r(t) + r* r(t) - r * r * r(t) H , (10) 

where a star denotes the convolution operation. We note that in this case 
the resolvent kernel is directly proportional to acceleration, so that r = 
and, by equation (Jjl]), k = for t < or r < r , i.e. before the acceleration 
is turned on. However, the character of memory that is indicated by k, 



F(r)=f(r)+ I" k(r-T')f(r')dr' 

J T 

PT — TO 

= /(r)+ / k(t)f(r-t)dt 



11 
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is more complicated than in case (1) due to the intricate relationship between 
r(t) and k(t) in equation fllPl). Even if the acceleration is turned off at r = n, 
it turns out that k does not vanish in general for r > T\ and could even be 
divergent; in fact, proving the latter point is the main purpose of this paper. 

Imagine, for instance, that k(t) is finite everywhere and decays exponen- 
tially to zero for t — > oo. Then in equation ([TT]) , as r — > oo long after the 
acceleration has been turned off at r = T\, the contribution of the nonlocal 
term in flTT| ) rapidly approaches a constant and we essentially recover the 
"nonpersistent" kinetic memory familiar from case (1). It turns out, how- 
ever, that in general case (2) involves situations with persistent or dynamic 
memory such that under certain conditions k(t) could diverge resulting in an 
asymptotically divergent F(t). 

The convolution (Faltung) type kernel is generally employed in many 
branches of physics and mathematics. As in equation (pT|) , to produce the 
nonlocal part of the output F(t), an input signal f(r — t) is linearly folded, 
starting from r and going backwards in proper time until To, with a weight 
function k(t) that is the impulse response of the system. The use of con- 
volution type kernels is standard practice in phenomenological treatments 
of the electrodynamics of media PTI, 12^, E3, feedback control systems p3 



etc. We find, however, that for the pure vacuum case the convolution kernel 
due to nonuniform acceleration in general leads to instability and is there- 
fore unacceptable. This proposition is proved in the following section for the 
translational and rotational accelerations of the observer. 

The simplicity of the kinetic memory versus dynamic memory has been 
particularly stressed by Hehl and Obukhov in their investigations of nonlocal 
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electrodynamics ^ZJ^ moreover, their work has led to the question of the 
ultimate physical significance of the convolution type kernel in the nonlocal 
theory of accelerated systems [ER E6[ . This question is settled in the present 



paper in favor of the kinetic memory kernel. 



4 Dynamic memory of accelerated motion 
4.1 Linear acceleration 

Imagine an observer at rest on the z- axis for —oo<t<tq. At r = tq, the 
observer accelerates along the positive ^-direction with acceleration g{r) > 0. 
For r > r , we set 



9{t) 



g(r') dr>, 



(12) 



TO 



C = cosh# and S = sinh^. The natural nonrotating orthonormal tetrad 
frame of the observer along its worldline is given by 



>) 

A (2) 



(C,0,0,S), Af 1} = (0,1,0,0), 
(0,0,1,0), Af 3) = (S,0,0,C). 



(13) 



In this case A(r) is given by 



A 



U V 
-V u 



u 



C 0' 
C 
1 



V = SI- 



S' 



(14) 



where Jj, {Ij)jk — ~^ijki is a 3 x 3 matrix proportional to the operator of 
infinitesimal rotations about the x*-axis. 

Let us first consider case (1), for which the kernel can be easily computed 
using equation (^), 



ko(r) = -g(r) 



h 
-h 



(15) 



so that when the acceleration is turned off at r = T\ the kernel ko vanishes 
with the acceleration for r > t%. On the other hand, k is simply constant 
for uniform acceleration (i.e. hyperbolic motion) with g(r) = go for r > 
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r . In the rest of this section, we focus attention on case (2) involving the 
convolution kernel. 

For the convolution kernel, the resolvent kernel is given, via equation (|9|), 

by 



r{r - r ) = g{r) 



SJ 3 Ch 
-Ch SJ 3 



(16) 



where (Jk)ij = $ij — Sikdjk- In principle, the convolution kernel can be com- 
puted via the substitution of equation (|Tj|) in equation fllOD ; however, this 
turns out to be a daunting task in practice. Imagine, for instance, that the 
acceleration is turned off at r = Ti, so that the resolvent kernel (|16"D has 
compact support over a time interval of length a = T\ — To and vanishes oth- 
erwise. It then follows that the r* n term in the expansion ( |T0| ) has compact 
support over a time interval of length na. The summation of series fllPf) turns 
out to be rather complicated, except for the case of uniform acceleration, i.e. 
g(r) = go for r > r , and the result is 



-go 



h 
-h 



(17) 



It is interesting to note that equation (|17D is the same as the result of case (1), 
equation (Jl5|), for uniform acceleration. 

In view of the difficulty of summing the series (|l0|) directly, we find it 
advantageous to use Laplace transforms, which we denote by an overbar, i.e. 
£{k(t)} = k(s), where 



k(s) 



e~ st k{t) dt- 



(18) 



then, taking the Laplace transform of equation ( |T0| ) and using the convolution 
(Faltung) theorem repeatedly, we arrive at 



k(s) = [I + f(s)]- 1 - I, 
which is consistent with the reciprocity between k and r. 



(19) 



4.2 Stepwise acceleration 

Let us specialize to a simple case of stepwise uniform acceleration, namely, 
we let g{r) = g for r < r < T\ and zero otherwise (see Figure |l|). In this 
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t 

9 



i 

1 1 

to n r > 

Figure 1: The linear acceleration of an observer that undergoes uniform 
acceleration go during a period a = T\ — t 2 of its proper time. If the area 
under the graph exceeds a critical value given by fi ~ 1.2931, then the 
convolution kernel leads to divergences. 



case, 



r(s) 



r\ r 2 
-f 2 fx 



(20) 



where fx(s) = q(s)J3 and f 2 (s) = p(s)Is. Here p(s) = £{gC} and q(s) = 
£{gS}. All 6x6 matrices that we consider in this paper have the general form 
(|20|), i.e. each is completely determined by two 3x3 matrices just as fx and 
fi characterize f in equation (|20|) ; we therefore write f — > [fx] f 2 ] to express 
this decomposition as in equation (p0|). To find the Laplace transforms of 
gC and gS, we note that in equation (|12"|), # = g (r — r ) for r < Tx and 
9 = Po = 9o( r i — To) f° r r > Tij therefore, 



p(s) ±g(s) 



9o 



sT9o 



[1 



-(s=F9o)"l 



(21) 



where a = ri — r = flo/go is the acceleration time interval. Using the results 
of Appendix B, we find from equation ( |I9"D that k(s) can be expressed as 



k(s)^[p Q(s)J 3 ]P P(s)I 3 
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where 



P{s) = -pM-e w + cosh/5 ) + A, sinh A)], 
= 71 [e tt (w sinh /3 - /3 cosh /3 ) + /3 ] . 



(22) 



(23) 



Here if := sa and the denominator D can be factorized as 



D = (we 



w 



(24) 



It is useful to recall that the kernel k — > [k\] k 2 ] refers to a system at rest on 
the z-axis for r < r that is uniformly accelerated at r = r with acceleration 
g until r + a = r 1; and then continues with uniform speed tanh/? along 
the positive z-direction for r > T\. Under certain conditions, it is possible 
to obtain series representations for k\ and k 2 (see Appendix C); however, to 
gain insight into the asymptotic behavior of k\ and k 2 it proves more fruitful 
to proceed with an investigation of the singularities of k\(s) = PoQ(s)Js and 
^2(5) = floP{s)h in the complex s-plane. This is due to a simple property 
of the Laplace transformation in equation (|1|) extended to the complex s- 
plane: let us suppose that the convolution kernel k(t) is a bounded function 
for all t = t — t > as one naturally expects of a function that represents 
memory; then, for any s in the complex plane with positive real part, i.e. 
Re(s) > 0, equation ( pT8[) implies that the absolute magnitude of k(s) should 
be finite, i.e. k(s) cannot be singular. Therefore, if we could show that k(s) 
has in fact pole singularities at complex values of s with Re(s) > 0, then it 
would simply follow that k(t) cannot be bounded for all t > and would 
thus be unsuitable to represent the memory of finite accelerated motion. 

We will prove the following result: If j3oexp(j3 ) > 3n/2, then the corre- 
sponding function k is unbounded for t > 0. It suffices to show that k has 
a pole in the right half of the complex s-plane. In fact, let us suppose that 
k has a pole at s = s , where Re(s ) > 0, but \\k\\ := sup t>0 \k(t)\ < 00. In 
this case, k has a pole in the half-plane 7i consisting of all complex numbers 
s such that Re(s) > | Re(so), and therefore \k\ is not bounded on H. On the 
other hand, for s e 7i, we have that 



in contradiction. Thus the rest of this subsection is devoted to the determi- 
nation of the poles of k(s) in the right half-plane. 
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The poles of k are elements of the zero set of D with w = sat and a > 0. 
Note, however, that the (real) zeros w = ±/3 are removable singularities. 
Poles in the right half-plane are the zeros of D with nonzero imaginary parts. 
In view of the definition of D, let us consider the complex roots in the right 
half-plane of the equation wexp(w) = b, where b is one of the real numbers 
±(3 exp(±j3 ). Because the zero set of this relation is symmetric with respect 
to the real axis, it suffices to consider only roots in the first quadrant of the 
complex w-plane. 

We set w = £ + ii], where £ > and 77 > are real variables, and note 
that wexp(w) = b if and only if 

= 5 cos 77, r\S = —5 sin 77. 

If this system of equations has a solution, then, by squaring, adding and 
rearranging, we have that rf = 6 2 exp(— 2£) — £ 2 or, since 77 > 0, 77 = 

V& 2 exp(-20-e 

There are several cases. For example, for b > 0, there is a pole in the 
right half-plane if the system of equations 

77 = \/b 2 e~ 2 t — £ 2 , = bcosrj 

has a solution with £ > and 77 mod 2n G (37r/2,27r). Similarly, for b < 0, 
there is a pole in the right half-plane if the system of equations 

77 = \/b 2 e~ 2 ^ — £ 2 , = bcosrj 

has a solution with £ > and 77 mod 27r e (77/2,77). 

A necessary condition for the relation 77 = a/6 2 exp(— 2£) — £ 2 to have a 
solution (£,77) is that £exp(£) < For b < 0, we must have £exp(£) < 
/3o exp(— /3 ); hence, there is a unique real number £ such that the necessary 
condition is met whenever £ < £ . On the other hand, for 6 > 0, the necessary 
condition, £ exp(£) < /3oexp(/3 ), is met if and only if £ < £ = A)- 

Let us view 77 as a function of £ and note that 7/(0) = 77(^9) = 0, and 

77^ = -e^6 2 - £ < 

for £ > 0. In particular, 77 decreases monotonically for < £ < £ - 

Consider the relation £exp(£) = bcosrj. At £ = 0, we have cos 77 = 0; 
therefore, the implicitly defined function 77 is such that 77(0) is an odd integer 
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Figure 2: The real branches of £ exp(£) = 5 cos 77 for 6 > 0. 

multiple of tt/2. At £o, we have cos// = ±1 according to the sign of b. In 
fact, 77 (£o) is an even multiple of n for 6 > and an odd multiple of it for 
6 < 0. Also, let us note that 

(£ + l)e^ = — b-j- sin 77. 

Suppose that b > 0. We will determine the positions of the real branches 
of the curve defined by £ exp(£) = bcosr/. For < 77 < tt/2, we have sin 77 > 
and drj/dl; < 0, so there is a real branch connecting the points (0, 7r/2) and 
(£0, 0) in the (£, r?)-plane. For n/2 < 77 < 3tc/2, we have cos r\ < 0; thus, there 
is no real branch in this region. There is a real branch connecting (0, 37r/2) 
and (£0, 27r) with drj/d^ > 0. This pattern continues as depicted in Figure 0. 
Note, however, that only the "increasing" branches correspond to poles in 
the right half-plane. Indeed, for b > 0, it is necessary that rj mod 2ir be in 
the interval (37r/2, 27r). In particular, the "lowest" branch corresponding to a 
pole connects the points (0, 3n/2) and (£0, 2tt). It is now clear that the curve 
defined by rj = \Jb 2 exp(— 2£) — £ 2 intersects an increasing branch with £ > 
if and only if b > 3n/2. The number of poles in the right half-plane increases 
by one as b increases past an odd multiple of 7r/2. 

Suppose that b < 0. In this case, the real branches of £ exp(£) = 6 cos 77 
exist only if cos rj < as in Figure |3] and a corresponding pole in the open 
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Figure 3: The real branches of £ exp(£) = b cos 77 for 6 < 0. 



right half-plane does not exist unless |6| > ir/2. Using the definition of b, this 
condition is equivalent to the requirement that /5 exp(— (3q) > n/2. But, the 
maximum value of /3oexp(— f3 ) is 1/e < 7r/2. Hence, negative values of 6 do 
not correspond to poles in the right half-plane. 

We conclude that the dynamic memory kernel k for stepwise uniform 
linear acceleration is unbounded for f3 = g a > 1.3. 



4.3 Rotation 

Imagine next an observer that is initially moving uniformly with speed v in 
the (x, y)-plane along a line parallel to the y-axis at x = po- At t = 0, x = po 
and y — 0, the observer starts rotating on a circle of radius po with uniform 
frequency Qo = v/po in the positive sense around the z-axis. Though the 
motion is continuous, there is no acceleration for t < and uniform circular 
acceleration for t > 0. The natural orthonormal tetrad frame of the uniformly 
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Figure 4: Schematic plot of the motion of the observer that undergoes 
stepwise uniform rotation of frequency O during a period a = T\ — r 2 of its 
proper time such that <p = jft a. If (po exceeds ir/2, then the convolution 
kernel leads to divergences. 



rotating observer is given by 

A/L = 7(1, — v sin <f, v cos ip, 0), 



(0, cos (p, sin ip, 0), 



A (i) 

A(2) = 7(1;, — sinyj, cos (p, 0), 



(25) 



A (3) 



(0,0,0,1) 



where 7 = (1 



"2 is the Lorentz factor and 9? = f2 t = 7^o r ; so that we 



have set r = in this case. Computing (j) a ^ for the tetrad frame (p5|), we find 
as expected that the translational acceleration has only a radial component 
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9i 



-vy 2 Qo and the rotational frequency is along the z-direction with 



magnitude f2 3 = 7 2 fi . Thus / = Af, where A — > [A x ; A 2 ] is given by 



Aa 



7 cosy? 7siny9 
— sin ip cos ip 
7 



A 5 



1 



— cos ip — sin <p 



(26) 



Let us first consider case (1); the kinetic memory kernel ko can be easily 
computed using the fact that for A given by equation fl26|) we have A -1 — > 
[Af;Ag]. Then we find that k -> [O J; -g • I], where fi= (0,0,7 2 fi ) and 
g = (— v 7 2 ^0; 0, 0) with respect to the orthonormal tetrad frame (E5T) . Thus 
fco is a constant kernel so long as the observer rotates uniformly; for instance, 
if the acceleration is turned off at T\ = a corresponding to ipo = 7^o«, then 
the observer will have uniform linear motion again with speed v for r > t\ 
and the kernel ko will vanish (see Figure ||). 

Let us now consider case (2); the dynamic memory kernel is given by the 
series (0) in terms of the resolvent kernel. This is given by equation @, 
r [r\; r 2 ], where 



r 1 = 7^0 



—7 sin (p 7 cos (p 
—7 cos ip — sin <p 
v 2 ^ 2 sm(p 



r 2 = t>7 2 fi 



7 sin ip 

cos ip 

7 sin (p — cos (p 



(27) 



The explicit calculation of k using the series fllP]) for the general case of step- 
wise uniform rotation from r = to t% = a is rather complicated; however, 
for n — >• 00 the calculation can be carried through and the result is a constant 
kernel given by k — > [O-I; — g ■ I]. Just as in the case of uniform translational 
acceleration (cf. Section 4), we have ko = k for uniform rotation as well. 

To calculate k for the stepwise uniform rotation of duration T 1 — T = a> 
0, we use Laplace transforms as in the previous section (see Figure |j). Let 
C = a~ 1 £{cos(y9} and 5" = a^^Cjsinyj}; then, with w = sa we find 



C ± iS' 



_ g-(i«T«¥>o) 



(28) 



19 



and hence the Laplace transform of the resolvent kernel is given by r 
[fx ; f 2 ] , where 



n = (f 



-j 2 S' 7 C" 
-jC -S' 
v 2 1 2 S' 



, r 2 = w7v?o 



' 1 S r 

o o a 

jS' -c o 



(29) 



Using methods given in Appendix B, equation (^) leads to 

k(s) - [ki(s)Ms)}, 

where 



7 2 Q lV 
-iP Q 
-v 2 ^ 2 Q 

- 7 Q" 
OOP 
7Q -V 



(30) 
(31) 



Here P and Q are given by 



V = -jj[w{-e w + cosy2 ) - V^osin^o], 
Q = — [e w (-wsm<p + v^ocos^o) - <A>], 



and the denominator P is given by 



£> = (we" - i<p e i(po ){we w + w"^ ). 



(32) 
(33) 

(34) 



It is interesting to note that if we formally substitute /3q for itpo in equa- 
tions (|32|)-(|3"4"D, we obtain results familiar from the previous subsection; 
specifically, under i<p — > j3 , V — > P, Q ^ iQ and P — > P, where P, Q 
and P are given in equations (|22])-(p4|). Therefore, the main results of the 
previous subsection can also be used in the analysis of stepwise uniform ro- 
tation; for instance, with appropriate modifications the explicit expressions 
given in Appendix C for the convolution kernel in a special case can be em- 
ployed here as well. However, since (po > 0, the singularities of P and Q are 
in general different from those in the previous subsection. 
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To determine the pole singularities of k(s) in the right half-plane in the 
case of stepwise rotation it suffices to consider the equation 

we w = i<p e iipo (35) 

with ip > 0. Indeed, note that if w is a solution of this equation, then the 
complex conjugate of w is a solution of wexp(w) = —iip exp(—i(p ). 

As before, let us set w = £ + %r\ and note that equation (|35|) is equivalent 
to the system of real equations given by 

iS- = Lp sm(i] - (f ), r/e e = ip cos(r? - (p ), (36) 

where £ > and <p > 0. We recall here that the solution w = i(p , i.e. 
£ = and i] = ip , of equations (|35D and (|36|) corresponds to a removable 
singularity. A necessary condition for system (^) to have a solution with 
£ > is that sin(r] — <p ) > and r]cos(r] — <p ) > 0; the latter condition 
means that cos(r/ — (po) and rj must have the same sign. 
Consider the system 

(£ 2 + r] 2 )e 2 ^ = <a> 2 , ^ = Vo sHv ~ fo). (37) 

If it has a solution (£,77), then it follows from display (|37| ) that 

rj 2 e 2li = ip 2 cos 2 (r] - <p ), 

and therefore 7/exp(£) = ±<p cos(r] — ip ). Comparing this result with sys- 
tem (|36|) , we conclude that we can use system (|37D for finding the poles if we 
keep in mind that rj and cos(^ — ifo) must have the same sign. 

The first equation in display ( |3~7D is equivalent to rf = ip 2 exp(— 2^) — £ 2 . 
Its graph in the right half-plane has the form depicted in Figure [5], where £0 
is the unique real solution of the equation £ exp(£) = (p Q . 

The poles we seek correspond to the intersections of the graph in Figure || 
with the real branches of the second curve in display ([37]). The intercepts 
of these branches with the 77-axis are given by the solutions of the equation 
sm(r] — ipo) = 0; that is, rj is equal to ip plus an integer multiple of ir. Along 
the line given by £ = £0, the intercepts are given by £0 exp(£ ) = (fo sm(rj—ipo). 
Because £oexp(£ ) = ¥0, these intercepts are the solutions of sm(r] — ip ) = 1; 
that is, rj is ipo + 7r/2 plus an integer multiple of 2ir. The shape of the 
branches connecting points on the two vertical lines (at £ = and £ = £0) 
is determined by the sign of cos(^ — <£> ) along the branch. Indeed, we have 
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Figure 5: The graph of rj 2 = ip 2 exp(— 2£) — £ 2 . 



already established that poles occur only at points where rj and cos(r/ — ip Q ) 
have the same sign. Note that 

(f + l)e ? = ip Q — cos(r? - <p ), 

and therefore the slope of the branch has the same sign as cos(r] — ipo). 
Moreover, only the branches with sin(?7 — tpo) > correspond to poles in the 
right half-plane. 

There are several cases depending on the size of ip . For < y? < 7r /2, it 
is easy to see that the important branches are as depicted in Figure These 
would not intersect the graph in Figure ^ hence, there are no poles in the 
right half-plane. 

We will next show that if <^ > 7r /2, then there is at least one pole in 
the right half-plane. For (po in this range, there is an integer j > 1 such 
that j7r/2 < tpo < (j + l)7r/2. In particular, we have that (po — J 7r /2 > 
and ipo — (j + l)7r/2 < 0. There are four cases. (1) Suppose that j is 
even and cos(j7r/2) = 1. The branch of the curve £ exp(£) = <^ sin(^ — <p ) 
with ^-intercept (po — jir/2 > has positive slope (like the upper branch 
in Figure [|). Because ip — jn/2 < ip , this branch intersects the curve 
depicted in Figure |5| in the upper half-plane. This point corresponds to a 
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Figure 6: The graph of £ exp(£) = <^ sin(r] — ip ) for < fo < tt/2. 



pole. Indeed, at the point of intersection sin(?7— (p ) > and rj cos(r/— ip ) > 0. 
(2) Suppose that j is even and cos(j7r/2) = — 1. The branch with //-intercept 
<Po — j'tt/2 > has negative slope and meets the line £ = £o with ordinate 

— (j + l) 7r /2 < 0. Hence, this branch intersects the curve depicted in 
Figure [| in the lower half-plane. This point corresponds to a pole. (3) 
Suppose that j is odd and cos((j + l)7r/2) = 1. The branch of the curve with 
77-intercept ipo — (j + l)vr/2 has positive slope and it meets the curve depicted 
in Figure |5| in the upper half-plane where the intersection point corresponds 
to a pole. For the subcase where j = 3 and <po = 3n/2, it is interesting 
to note that 77 = and £0, such that £oexp(£ ) = 37r/2, is the pole. (4) 
Suppose that j is odd and cos((j + l)vr/2) = —1. The curve with //-intercept 
^0 _ (j + l) 7r /2 has negative slope and —(po < tpo — (j + 1)^/2. Hence, this 
branch meets the curve depicted in Figure |5] in the lower half-plane where 
the intersection corresponds to a pole. 

We conclude that the dynamic memory kernel k for stepwise uniform 
rotation is unbounded for (p = ^VL^a > tt/2. 
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4.4 Smooth acceleration 



We have demonstrated that the convolution kernel k is unbounded for certain 
stepwise translational and rotational accelerations. Could this result be due 
to the discontinuities of these accelerations at tq and n? To prove that this 
is not the case, we are interested here instead in smooth accelerations that 
closely approximate the stepwise ones already studied. The translational 
and rotational cases are in fact closely related as we have demonstrated; 
therefore in this subsection we show the same result for the simpler case of 
smooth translational acceleration. 

Let us consider an acceleration g with compact support in the interval 
[tq, 7~i]. By the definition of A and the choice of g, the matrix A(r ) = A is 
the 6x6 identity matrix. Using this fact and equations @ and ([14]) , we find 
that r(t) = [g(r)S(r)J 3 ;g(r)C(r)I 3 ], where t = r-r , S(r) = sinhfl, C(r) = 
cosh^ and 9(t) is given by equation (|l^). It follows from equation (|18D that 
f(s) = [S(s)J 3 ;C(s)I 3 ], where 



C(s) 



POO 

±S(s) = e~ st g(t + T )e ±e{t+To) dt. 
Jo 



Using equation ( |T9| ) and the results of Appendix B, we find that the Laplace 
transform of the convolution kernel k is given by k(s) = [H\(s) J 3 ; 0,2(3)13], 
where 

has) = - l t s - -1, n 2 (s) 



[1 + S) 2 -C 2 ' zv ' (1 +S) 2 -C 2 ' 

We are interested in the zeros of the denominator 

(1 + S) 2 - C 2 = (1 + S + C)(l + S - C). 

It suffices to demonstrate that 1 + S + C has a zero in the right half of the 
complex s-plane. Because g has compact support in the interval [t ,ti], the 
function 1 + S + C is given by 



pa 

1+ / e- st e^ 9{a+To)da g(t + T Q )dt, 
Jo 



where a = r x — r . If (7 is the stepwise uniform acceleration considered 
previously, then this function reduces to 



^^(se as - g e a9 °: 

s- go 
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and, by the results in Subsection ^|for equation (|2~3|), if (3 exp((3 ) > 37r/2, 
it has a zero in the right half of the complex s-plane corresponding to a pole 
of k. We will show that such a pole persists for a smooth acceleration that 
is sufficiently close to the stepwise acceleration. 

For an arbitrary acceleration g with support in the interval [to,ti], we 
define the associated real-valued function ( on the interval [0, a] given by 
C(t) = g(t + r o)- Also, recall that the L 1 -norm of a real-valued function v 
defined on the interval [0, a] is given by 



v\\i : = / \v(t)\ dt. 



o 



Suppose that ( and v are real-valued functions defined on the interval 
[0, a] such that < oo and ||t>||i < oo, and consider the complex- valued 
analytic functions Z and T of the complex variable s given by 



pa 

Z(s) = 1+ / e- st eti c ^ dtT ({t)dt, 
Jo 

pa 

T(s) = 1+ / e- st eti v{a)da v{t)dt. 
Jo 



We will prove the following proposition. // Z has a zero in the open right- 
half of the complex s-plane and \\v — £||i is sufficiently small, then T has 
a zero in the open right-half of the complex s-plane. By a standard result 
from mathematical analysis (see |27]), if C is an L 1 function (for example, 



if ((t) = g(t + r ) for the stepwise acceleration g), then \\( — v\\i can be 
made as small as desired for a C°° function v. Hence, by the proposition, 
there is a smooth acceleration with compact support such that its associated 
convolution kernel is unbounded. 

Our proof begins with two estimates. For notational convenience, let us 
define 

((t) = efiCW^Cit), v(t) = eti v{u)da v(t). 
The first estimate is 

\T{s)-Z{s)\ < \\v-C\\i (38) 
for all s such that Re(s) > 0. To prove it, note that 

pa 

\T(s)-Z(s)\< / \e- st \\v(t) -£(t)\dt. 
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Because | exp(— st) \ < 1 for Re(s) > 0, we have the inequality 

|T( a )-Z( S )|<||t)-C||i 
for all s in the closed right half-plane. The second estimate is 

||0_{|| 1 < e IICIIi( 1 + a || C ||) e ll«-Clli|| t; _ C || 1 . (39) 

To prove it, we have the triangle law estimate 

\v(t)-C(t)\ < \eti v ^ da v{t) - eti v{cT)da ({t)\ 

+ \eti v{a)da C{t) ~ e /( ' c(ff)a!a CWI 
< e /oK <J )l (iCT | u - (| + \( > \\e^ v ^ )da - e ti c - {a)da \ (40) 

and, by the mean value theorem (applied to the exponential function), the 
inequality 



Jo Jo 

where q is some number between f*v(a)da and LQ{a)da. If <r < 0, then 
exp(^) < 1; and if <; > 0, then q < max{||t>||i, ||C||i}- Hence, 

e ? < e max {IMIi,IICIIi} 

and, because \\v\\i < \\v — C||i + HCIIi? we have that 

e < < e llC||i e l|t;-CI|i. 

Using this result and the estimate (0), it follows that 

\v(t) - C(t)\ < e^\v(t) - C(t)\ + HCHel^ 1 ^ 1 J" \v{a) - C(a)\ da 
< ell^ e ll^(i,;(t)-C(t)i + ||Clllb-Cl|i)- 

Therefore, 



|t> — CIU = / \v~C\dt 
Jo 

< e IICBi e ll«-flli(|| t; _ C || 1 + a || C |||| t; _ C || 1 ; 
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and a rearrangement of the right-hand side of the last inequality gives the 
desired result. 

In the rest of this section, we let ( represent the stepwise uniform linear 
acceleration and v the smooth linear acceleration that approximates it suffi- 
ciently closely. We then choose a circle centered at a zero of Z in the open 
right half of the complex s-plane such that the circle does not pass through 
a zero of Z and such that the circle is contained in the open right half-plane. 
Let k, a complex- valued function defined on the interval [0, 2n], be a contin- 
uous parametrization of this circle and define two new functions Kz and /tx 
on this interval by 

k z (0) = k t (0) = T(/c(0)). 

The images of these functions are closed curves in the complex s-plane. In 
complex analysis, the principle of the argument theorem ]25] for an analytic 



function A relates the winding number of the image of ka with respect to 
the origin to the number of zeros of the function A inside the circle, provided 
that the circle does not pass through any zero of A. If we show that Kz and 
kt are homotopic and therefore have the same winding number with respect 
to the origin and that the circle does not pass through a zero of T, then Z 
and T must have the same number of zeros inside the circle. 

We claim that if ||i> — Clli is sufficiently small, then the image of k does 
not pass through a zero of T. To prove the claim, note that 

m := mm{\K Z ($)\ : < <& < 2vr} > 

(because k does not pass through a zero of Z) and using the triangle inequal- 
ity 

< m < |«t(i?)| + || «t — K z\\, 
where ||kt — Kz\\ is the supremum of |/ct($) — K z{$)\ for < $ < 2tt. Using 



the estimates (|38D and (0), we have that 



M0) - «zWI < e^K^(l +a\\(\\)\\v - C||x. (41) 

By the estimate \\Ky — k z \\ can be made small, say less than m/2, by 
taking \\v — C||i sufficiently small. For all v satisfying this requirement, which 
we impose for the remainder of the proof, we have that |ky($)| > 0; that is, 
k does not pass through a zero of T. 
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It remains to show that is homotopic to Kz- Assuming this homotopy 
relation, the image curves of and Kz would have the same winding number 
with respect to the origin. By the choice of k and the argument principle 
(see [0), the curve Kz has a nonzero winding number. Hence, kx would 
have the same nonzero winding number. Again, by the argument principle, 
T must then have a zero in the disk bounded by the circle parametrized by 
k, which is the desired result. 

To complete the proof we need to show that kt and kz are indeed 
homotopic. Let C denote the complex numbers. We will show that H : 
[0, 1] x [0, 2vr] -> C \ {0} given by 

#(a,i?) = k z (#) - ct(k z (&) - kt(0)) 

is the required homotopy. By inspection, H is continuous, if(0,i?) = k.z(i?) 
and if = Ky^). Hence, it suffices to show that if(<r, i?) ^ for all 
(a, i?) G [0, 1] x [0,27r]. By our choice of v, we have that ||k t — k z \\ < m/2; 
therefore 

\H(a,0)\ > |« z (i?)| - \<t\\kt(0) - kz(&)\ >m- ||« T - «z|| > m/2, 
as required. 

We conclude that the dynamic memory kernel k for the smooth linear 
acceleration that closely approximates the stepwise acceleration is unbounded 
if the area under the graph of g(r) exceeds a critical value ~ 1. 

5 Discussion 

We have investigated the properties of the nonlocal kernel that is induced by 
accelerated motion in Minkowski spacetime. The physical principles outlined 
in this paper do not completely determine the kernel; therefore, simplifying 
mathematical assumptions need to be introduced in order to identify a unique 
kernel. Two possibilities have been explored in this work corresponding to 
kinetic memory (ko) and dynamic memory (k). We show that for accelerated 
motion that is uniform (linear or circular), the two kernels give the same con- 
stant result ko = k. They differ, however, if the acceleration is turned off at 
a certain moment. We have therefore studied piecewise uniform acceleration 
(linear and circular) and have demonstrated that the dynamic memory (con- 
volution) kernel could be divergent and is therefore ruled out. Furthermore, 



28 



this conclusion is shown to be independent of the stepwise character of the 
linear acceleration considered. 

The use of convolution kernels is standard practice in the nonlocal electro- 
dynamics of continuous media, where it is assumed phenomenologically that 
memory always fades. In our treatment of acceleration-induced nonlocality 
in vacuum, however, the behavior of memory must be determined from first 
principles. In this connection, the possible advantage of kinetic memory in 
terms of its simplicity was first emphasized by Hehl and Obukhov p5, 26fl . 

The theory developed here is applicable to any basic field; however, for 
the sake of concreteness and in view of possible observational consequences, 
we employ electromagnetic radiation fields throughout. A basic consequence 
of the nonlocal theory of accelerated systems is that it is incompatible with 
the existence of a basic scalar field; that is, in this case A(r) = 1, ko — 
and the nonlocality disappears so that a basic scalar radiation field can stay 
completely at rest with respect to a rotating observer in contradiction with 
our fundamental physical assumption. This prediction of the nonlocal theory 
is in agreement with present experimental data. Further confrontation of the 
nonlocal theory with observation is urgently needed. 



Appendix A 

Consider an integral equation of the form 

<f>{x) = il){x) + e I K{x,y)<j){y)dy, (Al) 

J a 

where ip is a continuous function, the kernel K is continuous and e is a 
constant parameter. There is a unique continuous resolvent kernel R such 
that 

px 

r()(x) = <p(x) + e R(x,y)iP(y)dy. (A2) 

J a 

In turn, K can be thought of as the resolvent kernel for R; this follows from 
the complete reciprocity between K and R. 

The proof of the existence and uniqueness of the resolvent kernel is by 
successive approximation. In fact, the solution <p can be obtained as the 
uniform limit of the sequence of continuous functions {4> n }'^Lo defined as 
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follows: (f>o(x) = if)(x) and 

PX 

(pn+i(x) = ip(x) +e K(x,y)(f) n (y) dy. (A3) 

J a 

Thus 

PX 

4>i{x) = tf;(x) + e / K{x,y)^{y)dy, (A4) 

J a 

px py 

(f) 2 (x) = tp(x) + e / K(x,y){ip(y) + e I K(y, z)ip{z) dz] dy 

(A5) 

The integration in ( |A5|) is over a triangular domain in the (y, z)-plane defined 
by the vertices (a, a), (x, a) and (x, x). Changing the order of the integration 
in (|A5|) results in the equality 



(A6) 



<pi{x) + e 2 I K(x,y) K(y,z)ip(z)dzdy. 



K(x,y) I K(y,z)ip(z)dz dy 

' a 

K(x,y)K(y,z)dy ip(z) dz. 
Let us define the successive iterated kernels of K by Ki(x, z) = K(x, z) and 

K n+ x(x,z) = J K(x,y)K n (y,z)dy. 
Then we can write ( |A5| ) as 

PX 

4>2{x) = <j>i(x) + e 2 / K 2 (x, z)ip(z) dz, 

J a 

and similarly 

PX 

fai?) = fa(x) + e 3 / K 3 (x, z)tp(z) dz, 

J a 

etc., such that in general 

PX 

(j>m(x) = (j) m -x{x) + e m K m (x, z)^(z) dz. (A10) 



(A7) 



(A8) 



(A9) 



30 



Iterating ( |A10| ) for m = 1, 2, 3, . . . , n and summing the equations results 



in 



(f) n (x) = lf)(x) + 

which can be rewritten as 

i[)(x) = 4> n {x) + e 



J2e m K m (x,z) 



m=l 



ip(z) dz, 



(All) 



m=l 



4>{y) dy. (A12) 



It can be shown that the uniform limit as n — ► oo exists (see [|10], [TT], |T^] 
Thus, we obtain equation ( |A"2j ) with 



(A13) 



R(x,y) = -Y,t n ~ 1 K n (x,y). 



n=l 



In case (1), K(x,y) = k (y), the iterated kernels K n for n > 1 and the 
resolvent kernel i? are in general functions of both x and y. 

In case (2), K(x,y) = k(x — y), i.e. the kernel is of the convolution 
(Faltung) type, it follows from ( |A7|) that 



(A14) 



where x — z = t and x — y = u; therefore, all of the iterated kernels are of the 
convolution type and can be obtained by successive convolutions of k with 
itself. More precisely, let a star denote the Faltung operation, 



k n+ i(t) = / k(u)k n (t - u) du, 



<l>*x{t)= (f>{t)x{t-u)du = x*<f>{t), 



(A15) 



and write 0* 2 = <p*4>, etc. Then, the resolvent kernel ( |A13| ) can be expressed 
as R(x, y) = r(x — y), where 



it) = -^V- 1 *;*"^). 



(A16) 



n=l 
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Appendix B 



In this paper, we deal with 6x6 matrices of the form 

M = 



A B 
-B A 



(Bl) 



where det A ^ and det B = 0. The inverse of the matrix M. is given by 



M 



-i 



G H 
-H G 



where 



G= (A + BA- L B)-\ H = —GBA 



-l 



A~ l BG. 



fB2l 



(B3) 



Appendix C 

Let us rewrite P(s) and Q(s) given by equations and (|23"1) in the form 
2(3o P(s) = f- - — -sl, (CI) 



2/3o Q(s) = -2 + f- + — (C2) 



i - C+ i + c 

where w = sa, /3 = go a an d C± are given by 

A) 







exp(-w ± Pq) 



(C3) 



If we assume that Re(s) > go> then \(±\ < 1. We can therefore expand 
(1 =F C±) _1 in powers of C± and use the relation 



£ \ u na (t) 



(t — na) 
('-!)! 



e-i 



(C4) 



for integers n > and £ > 1 to find fc(t) — > [fci(t); ^(i)]- Here we use unit 
step functions such that u na {t) = u (t — na) and u (t) is the standard unit 
step function, i.e. u (t) = 1 for t > and u (t) = for t < 0. 
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We find that k\ = &4J3, k 2 = k 2 I 3 and 

00 ^1(0 = SlU a (t) + C 2 U 2a (t) + S 3 U 3a (t) + C 4 M 4 a(t) 
Q^fait) + «o(0 = Cili a (*) + S 2 U 2a (t) + C 3 U 3a (t) + S 4 M 4 a(t) 

where 



Cn i 5 n 



(got - n(3 ) n - 1 _ (g t - n(3 ) n 



(n-1)! 



(C5) 
(C6) 



(C7) 



Note that for any fixed value of t, only a finite number of terms contribute 
to the kernel k(t). 
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